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Abstract 



(-H Dirac structures unify both presymplectic and Poisson structures, which serve an 

4—^ implicit generahzation of Lagrangian and Hamiltonian systems including the case of 

non-holonomic mechanics as well as the case of degenerate Lagrangians. It is known 
that a Dirac structure represents a power conserving interconnection structure between 
' ' physical systems. In this paper, we investigate the interconnection of distinct Dirac 

^ structures and associated physical systems. First, we make a brief review on induced 

^ Dirac structures and Lagrange-Dirac dynamical systems. Second, we consider how 

ly-^ distinct Dirac structures D\,...,Dn can be interconnected through a Dirac structure 

-Dint- To do this, we introduce a tensor product called the bowtie product, [xi, of 
7-H Dirac structures and then show that the interconnection of Dirac structures can be 

given by {Di © • • • © _D„) cxi Dint. We also explore variational structures associated to 
jy-^ the interconnection of Lagrange-Dirac systems. Lastly, we demonstrate the theory of 

interconnection of Dirac structures and associated Lagrange-Dirac dynamical systems 
by some examples including electric circuits, nonholonomic mechanical systems, and 
( simple mass-spring mechanical systems. 



^ Contents 

1 Introduction 2 

2 Dirac Structures 4 

3 Lagrange Dirac Dynamical Systems 6 

4 Interconnection of Dirac Structures 9 

5 Interconnected Lagrange-Dirac Dynamical Systems 17 



* Research partially supported by JSPS Grant-in-Aid 20560229, JST-CREST and Waseda University 
Grant for SR 2010A-606. 



1 



1 Introduction 2 



6 Examples 22 

7 Conclusions and Future Works 29 
A Appendix: The Bowtie Product [xi and the Dirac Tensor Product Kl 29 



1 Introduction 

A large class of physical and engineering problems can be described in terms of Lagrangian 
and Hamiltonian systems, such as controlled mechanical systems, externally forced and dis- 
sipative systems, collisions, electric circuits, stochastic systems, and field theories such as 
electromagnetism and elasticity. Such systems are usually categorized into Lagrangian or 
Hamiltonian systems with constraints, whether holonomic or nonholonomic; however, one 
may be faced with serious problems on ever larger scales becoming more and more heteroge- 
neous, involving a mixture of mechanical and electrical components, with flexible and rigid 
parts, and elastic or magnetic couplings (see, for instance, Yoshimura [1995]; Bloch [2003] 
and Afshari, Bhat, Hajimiri and Marsdcn [2006]). Dirac structures are known as a powerful 
tool for providing a natural geometric framework for describing the structure of such diverse 
Lagrangian and Hamiltonian systems on the product of configuration spaces by introduc- 
ing nontrivial constraints that model interactions between subsystems, which is called an 
interconnection. The notion of an interconnection was originally developed by Kron [1963] 
in his book of "Diakoptics" for modeling complicated physical systems. The word "diakop- 
tics" denotes a procedure of tearing a dynamical system into disconnected subsystems as 
well as interconnecting them to reconstruct the original system. The original system may 
be regarded as an interconnected system of subsystems, with constraints induced from the 
interconnection. Interconnection has been studied through the lens of power conservation 
motivated by the bond graph theory established by Payntcr [1961]. In electric circuit theory, 
the interconnection network has been widely employed through nonenergic multiports (see 
Brayton [1971]; Wyatt and Chua [1977]), where Kirchhoff's current law gives a constraint 
distribution on a configuration charge space such that Tellegen's theorem holds. In mechan- 
ics, it was shown by Yoshimura [1995] that kinematical constraints due to mechanical joints, 
coordinate transformations as well as force equilibrium conditions in dAlcmbcrt principle 
can be represented by the interconnection modeled by nonenergic multiports using the bond 
graph theory in the context of Lagrangian mechanics. 

It was demonstrated by van dcr Schaft and Maschkc [1995] that the interconnection can 
be represented by Dirac structures induced from Poisson structures and also that nonholo- 
nomic systems and L-C circuits can be represented by implicit Hamiltonian systems. On the 
Lagrangian side, it was shown by Yoshimura and Marsdcn [2006a] that nonholonomic me- 
chanical systems and L-C circuits (as degenerate Lagrangian systems) can be formulated by 
implicit Lagrangian systems associated with the Dirac structures on the cotangent bundles 
induced from given constraint distributions. From the view point of variational structures, 
it was shown by Yoshimura and Marsdcn [2006b] that the standard implicit Lagrangian 
systems, namely, implicit Euler- Lagrange equations for unconstrained systems can be de- 
rived from the Hamilton-Pontryagin principle and also that the Lagrange-Dirac dynamical 
systems can be formulated in the context of Lagrange-d'Alembert-Pontryagin principle. 

From the viewpoint of control theory, a notion of implicit port- controlled Hamiltonian 
(IPCH) systems (an implicit Hamiltonian system with external control inputs) was devel- 
oped by van dcr Schaft and Maschkc [1995] (see also Bloch and Crouch [1997], Blankcnstcin 
[2000] and van dcr Schaft [1996]) and much effort has been devoted to the passivity based 
control for the interconnected mechanical systems (Ortega ct al [1998]). For the case of 
regular Lagrangian systems, the equivalence between controlled Lagrangian (CL) systems 
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and controlled Hamiltonian ( CH) systems was shown by Chang, Bloch, Leonard, Marsden, 
and Woolsey [2002] and the reduction theory for the CL systems was shown by Chang and 
Marsden [2003]. For the general case in which a given Lagrangian is degenerate, an im- 
phcit Lagrangian analogue of IPCH systems, namely, a notion of implicit port- controlled 
Lagrangian (IPCL) systems for electrical circuits was constructed by Yoshimura and Mars- 
den [2006c] and Yoshimura and Marsden [2007a], where it was shown that L-C transmission 
lines can be represented in the context of the IPLC system by employing induced Dirac 
structures with Lagrange multipliers. 

Recently, it was shown by Ccrvera, van dcr Schaft, and Banos [2007] that any power- 
conserving interconnection of IPCH systems is another port-Hamiltonian system, with the 
Dirac structure being the composition of the Dirac structures of its constituent parts, Hamil- 
tonian the sum of the Hamiltonians, and resistive relations determined by the individual 
resistive relations. As a control application, the feedback interconnection of a plant port- 
Hamiltonian system with a controller port-Hamiltonian system can be represented by the 
composition of a plant Dirac structure with a controller Dirac structure. 

Goals and Organizations of the Paper. One of the main goals in this paper is to 
develop a simple expression for the interconnection of distinct Dirac structures and their 
associated mechanical systems; namely, we will show how to obtain an interconnection 
of distinct Dirac structures _Di,...,Z?„ on distinct manifolds Mi,...,M„ as another Dirac 
structure D on Mi x • • • x M„ by "interconnecting" _Di, D„. To do this, we will introduce 
two mathematical ingredients, namely, the Dirac sum and the howtie product ixi of Dirac 
structures and we will show how two distinct Dirac structures can be interconnected to yield 
another Dirac structure via an interconnection Dirac structure (usually labeled ZJjnt in this 
paper) under an assumption of transversal intersections. In particular, the interconnection 
of Dirac structures \s D = (Di® • • -(BDn) cxi Mnt- The interconnection Dirac structure Z3i„t 
provides an interface to interconnect "distinct" Dirac structures Di, D„ to yield another 
Dirac structure D. We will also show how the associated interconnection of Lagrange- 
Dirac dynamical systems can be developed in the context of variational structures. We will 
demonstrate our theory by some examples of circuits and mechanical systems. 

This paper consists of the following parts: 

1. We briefly review Dirac structures and we show how associated Lagrange-Dirac dy- 
namical systems can be developed by the Lagrange-d'Alembert-Pontryagin principle 
(see Yoshimma and Marsden [2006b]). 

2. We show how distinct Dirac structures can be interconnected to yield another Dirac 
structure on the product manifold through an interconnection Dirac structure. To do 
this, we develop a map called the Dirac sum and a tensor product called the bowtie 
product [XI between Dirac structures. 

3. We explore the variational structure of interconnection of Lagrange-Dirac dynamical 
systems; namely, the interconnected Lagrange-Dirac dynamical system can be formu- 
lated in the context of the Lagrange-d'Alembert -Pontryagin principle. 

4. We demonstrate some illustrative examples of circuits and nonholonomic systems, 
which can be generalized to more complex systems. 

Acknowledgements. We are also very grateful to Henrique Bursztyn, Mathieu Desbrun, 
Melvin Leok, Tomoki Ohsawa, Arjan van der Schaft and Alan Weinstein for useful remarks 
and suggestions. The research of H. J. is supported by Caltech and JST-CREST in part. 
The research of H. Y. is partially supported by JSPS Grant-in- Aid 20560229, JST-CREST 
and Waseda University Grant for SR 2010A-606. 
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2 Dirac Structures 

Linear Dirac Structures. First, let us recall the definition of a linear Dirac structure, 
namely, a Dirac structure on a vector space V; we assume that V is finite dimensional for 
simplicity (see, Courant and Wcinstein [1988]). Let V* be the dual space of V, and (• , •) be 
the natural paring between V* and V. Define the symmetric paring ((•, •)) on V (BV* by 

{{{v,a),{v,a))) = {a,v) + {a,v), 

for (w,a), (-0, a) &V ®V*. 

A Dirac structure on ^ is a subspace D C V (B V* such that D = , where is the 
orthogonal of D relative to the pairing ((•, •)). 

Dirac Structures on Manifolds. Let Af be a smooth manifold^ and let TM T*M 

denote the Whitney sum bundle over M, namely, the bundle over the base M and with 
fiber over the point x G M equal toT^M xT*M. A subbundle D C TM ®T*M is called 
an almost Dirac structure on M, when D{x) is a Dirac structure on the vector space T^M 
at each point x G M. A given two-form on M together with a distribution A on M 
determines an (almost) Dirac structure on M as follows: For each x € M, 

D{x) = {{v, a) e T^M X T*M \ v € A{x), and 

a{w) — nA{x){v,w) for all w E A{x)}, 

where Ha is the restriction of fl to A. We can also restate the Dirac structure as, for each 
X e M, 

D = {{v,a) e T^M X T*M \ v e A(x) and n\x)v - a e A°(x)}, 
where A° is an annihilator of A. 

We call D an integrable Dirac structure if the following integrability condition 

{£x,a2,X3) + (£x,a3,^i> + {fx.ai^X^} = (2.2) 

is satisfied for all pairs of vector fields and one-forms (Xi.ai), {X2,a2), (X^^as) that take 
values in D, where £x denotes the Lie derivative along the vector field X on M. 

Remark. Let T{TM®T*M) be a space of local sections of TM ®T*M, which is endowed 
with the skew-symmetric bracket [ , ] : T(TM ®T*M)x T{TM ®T*M)^ T{TM ® T* M) 
defined by 

[(Xi,ai), {X2,a2)] := {[Xi,X2] , £xia2 - £x2<^i + d {a2,Xi)) 
= ([-'^i, -'^2] , ixidaa - ixadai +d{a2,Xi)). 

This bracket is the one that Courant [1990] originally developed and it does not necessarily 
satisfy the Jacobi identity. It was shown by Dorfman [1993] that the integrability condition 
of the Dirac structure D C TM ® T* M given in equation (2.2) can be expressed as 

[r(i?),r(i?)]cr(i?), 

which is the closedness condition of the Courant bracket (see also Dalsmo and van der Scliaft 
[1998] and Jotz and Ratiu [2008]). 

The following Lemma is useful and the proof may be found in Yoshimura and Marsdcn 
[2006a] and Courant [1990]. 

Lemma 2.1. A subbundle D C TM © T* M is isotropic with respect to ((•, •)) if and only if 
{{{v,a),{v,a)))=0,y{v,a)eD. 

^In this paper, all the geometric objects such as manifolds, vector bundles and sections are assumed to 
be smooth. 
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Induced Dirac Structures. One of the most important and interesting Dirac structures 
in mechanics is an induced Dirac structure from kinematic constraints, whether holonomic 
or nonholonomic. Such constraints are generally given by a distribution on a configuration 
manifold (as to the details, refer to Yoshimura and Marsden [2006a]). 

Let Q be a configuration manifold. Let TQ and T*Q be the tangent and cotangent 
bundles of Q. Let Aq C TQ be a regular distribution on Q and define a lifted distribution 
on T*Q by 

At*q = (Tttq)-' (Aq) C TT*Q, 

where ttq : T* Q Q is the cotangent projection so that its tangent is a map Tttq : TT*Q — > 
TQ. Let n be the canonical two-form on T*Q. The induced Dirac structure D on T*Q is 
the subbundle of TT*Q T*T*Q, whose fiber is given for each {q,p) e T*Q as 

D{q,p) = {{v,a) G T^g,p)T*Q x Tl^^p^T*Q \ v e At.q(9,p), and 

a{w) = nAQ{q,p){v,w) for all w e AT*Q{q,p)}. 

This is, of course, a special instance of the construction (2.f). Notice that the induced Dirac 
structure can be restated by using the bundle map il^ : TT*Q T*T*Q as follows: 

D{q,p) = {{v,a) e T^,,p)T*Q x Tl^j,^^T*Q \ v e AT*Q(g,p), and 

a-n\q,p)-veA°T,Qiq,p)}, 

where A^-.g is the annihilator of At*q. 

Remark. If there exists no constraint, namely, the case in which Ag — TQ, the Dirac 
structure D corresponds to the canonical Dirac structure, which is given by a graph of 
the bundle map 51^ associated to Q. 

Local Expressions. Let be a model space for Q and U be the range in ^ of a chart 
on Q. Then, the range of the induced chart on TQ is represented hy U x V, while T*Q 
is locally given by J7 x y*. In this representation, {q,q) are local coordinates of TQ and 
{q,p) are local coordinates of T*Q. Further, TT*Q is represented by {U x V*) x (V^ x V*), 
while T*T*Q is locally given hy (U xV*) x (V* xV). In this local representation, iq,p, q,p) 
are the corresponding coordinates of TT*Q, while {q,p,/3,v) are the local coordinates of 
T*T*Q. 

Using ttq : T*Q Q locally denoted by {q,p) q and Tttq : {q,p,q,p) i-^ {q,q), it 
follows that 

At'Q = {w = {q,p, q,p)\q(zU,qe A{q)} 
and the annihilator of Ayg is locally represented as 

A?^.Q ^{a = {q,p, p,v)\q^U,P(, A° (q) and v = 0} . 

Since we have the local formula Q'^{q,p) ■ w = {q,p, —p, q), the condition a — {q,p) ■ w G 
A^-.g reads (3 + p e A°{q) and v — q — 0. Thus, the induced Dirac structure is locally 
represented by 

D{q,p) = {iiq,p),{P,v))\qeA{q),v^q, /3+peA°(g)}, (2.3) 
where A°{q) C T*Q is the polar of A{q) cTqQ. 
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3 Lagrange— Dirac Dynamical Systems 

In this section, let us recall the definition of implicit Lagrangian systems or Lagrange-Dirac 
dynamical systems by following Yoshimura and Marsdcn [2006a, b] . 



Lagrange— Dirac Dynamical Systems. Let L : TQ -> M be a Lagrangian, possibly 
degenerate. Let Dhe a Dirac structure on T*Q induced from a given distribution Aq C TQ. 

Recall that a partial vector field X : TQ (B T*Q TT*Q is defined as a mapping 
such that tt*q o X ~ pr^.g, where pr-^.g : TQ © T*Q — > T*Q; {q,v,p) i~7> i<l,p) and 
tt'q : TT*Q — > T*Q; {q,p,q,p) {q,p) is the tangent projection.^ The partial vector field 
X -.TQ® T*Q TT*Q is a map that assigns to each point (g, v,p) £TQ® T*Q, a vector 
in TT*Q at the point {q,p) G T*Q; we write X as 

X{q,v,p) = {q,p,q,p), 

where we note that q = dq/dt and p — dp/dt are understood to be functions of (g, 

Let El ■ TQ T'*(5 — > M be the generalized energy given by 

EL{q,v,p) := {p,v) - L{q,v). 

An implicit Lagrangian system or a Lagrange Dirac dynamical system is the 

triple {El,D, X) that satisfies, for each {q, v,p) G TQ®T*Q and with P = FL(Aq), namely, 
{q,p) = {q,dL/dv), 

{X{q,v,p),dEL{q,V,p)\T^^^^^p) e D{q,p), (3.1) 

where FL : TQ T*Q is the Legendre transformation. In the above, dEi^{q,v,p) : 
T{q,v,p){TQ®T*Q) ^- M is a linear function on T(q.v,p)iTQ ®T*Q) at each point {q,v,p) G 
TQ © T*Q and hence the restriction dE^iq, v,p)\T^g ^ {—dL/dq, v) may be regarded as 
a linear function on T^q p-^P such that {q,p) = {q, dL/dv) E P = F_L(Aq). 



Local Expressions. It follows from (2.3) that the Lagrange-Dirac dynamical system 
{X,dEL\Tp) E D is locally given by 

P^TT' ^^i-eAgq, p- — EA°Q{q). 
ov oq ^ 

Remark. If a partial vector field X{q,v,p) — {q,p,q,p) satisfies the condition for a 
Lagrange-Dirac dynamical system {X,dEL\Tp) E D, then the Legendre transformation 
iliP) = {Qt dL/dv) E P — FL(Aq) is consistent with the equality of base points and the 
condition itself gives the second-order condition q = v. In other words, the partial vector 
field of the Lagrange-Dirac dynamical system is uniquely given on the graph of the Legendre 
transformation. 

For the case in which no kinematic constraint is imposed, i.e., Ag — TQ, we can develop 
the standard implicit Lagrangian system, which is expressed in local coordinates as 

dL dL 

P = l^v, p = — , 

ov aq 

which we shall call implicit Euler Lagrange equations. Note that the implicit Euler- 
Lagrange equations include the Eulcr-Lagrangc equations p = dL/dq, the Legendre trans- 
formation p = dL/dv as well as the second-order condition q^ v. 



^Compare this to the definition of a vector field over T*Q as a mapping Y : T*Q — ^ TT*Q such that 
TT*Q oY = Idy.Q 
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Example: Harmonic Oscillators. We illustrate the implicit Euler-Lagrange equations 
by the example of a linear harmonic oscillator. In this case, the configuration space is 
Q = M representing the position of a particle on the line. The Lagrangian is given by 
L = v'^ /2 — /2. We find the generalized energy to be EL{q,v,p) = pv — + q'^/2. The 
canonical Dirac structure D{= graph(f2'')) can be easily constructed by using the canonical 
symplectic structure f2 on T*Q. The Lagrange-Dirac system is given by {El,D,X), which 
satisfies 

{{q,p),dELiq,v,p)\Tp) G D{q,p), 

and it follows dEL{q,v,p)\TP — i}^{q,p) ■ {q,p). We see dEL{q,v,p)\TP = vdp + qdq and 
f2''(g,p)(g,p) = —pdq + qdp. So, the dynamics is equivalent to the following equations: 

q^v, p=-q, 

together with the Legendre transform p = v. 

Lagrange-Dirac Dynamical Systems with External Forces. Here, we consider the 
case of a system with external force fields. Let prg : TQ © T*Q — Q; {q,v,p) n- q be 
a natural projection. An external force field F : TQ (B T*Q — )■ T*Q induces a horizontal 
one-form F onTQ®T*Q such that, for (g, v,p) e TQ ® T*Q, 

F{q,v,p) ■ w = {F{q,v,p),TwQ{w)) , 

where F = {q,v,p, F{q,v,p), 0,0) E T^^^^ (TQ ® r*Q) and w = {q,v,p,6q,Sv,Sp) e 
Ti<j,y,p){TQ(BT*Q). 

Given a Lagrangian L : TQ — > M (possibly degenerate), a Lagrange-Dirac dynamical 
system with an external force field is defined by a quadruple {El,D,X,F), which 
satisfies, in coordinates {q,v,p) £ TQ T*Q, 

{X{q,v,p)AdEL - F)iq,v,p)\Tp) e Diq,p). (3.2) 

It follows that the local expression of a Lagrange-Dirac system in equation (3.2) may be 
given by 

dL 

q = V, p-^-F{q,v,p)eA°Qiq), qeAQiq) (3.3) 

and with 

_ ^ 
^ dv 

The curve {q{t),v{t),p{t)), h < t < t2 in TQ ® T*Q that satisfies the condition (3.2) is 
a solution curve of {El,D, X, F). 

Power Balance Law. The time derivative of EL{q,v,p) ~ {p,v) — L{q,v) restricted to 
the solution curve {q(t),v{t),p{t)) reads 

^^EL{qit),v{t),p{t)) = {F{q{t),v{t),p{t)),qit)) , 
where we have employed p{t) = {dL / dv){t) and q{t) = v{t). 
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Example: The Harmonic Oscillator with Damping. As before, consider the har- 
monic oscillator in the setting of Q = M, L = — = pv — + and D{= 
graph(f]'')) and add a damping force F : TQ®T*Q T*Q defined by F{q, v,p) = -{rv)dq, 
where r G is a positive damping constant. Then, F = (q, rw, 0, 0). The formulas in 
equation (3.3) give us the equations: 

q = V, p + q + rv = 

with the Legendre transformation p = v. 



The Hamilton-Pontryagin Principle. Lagrange-Dirac dynamical systems can describe 
general scenarios such as non-holonomic degenerate Lagrangian systems because of a firm 
basis in variational structures. Given a Lagrangian L : TQ — > M, unconstrained equations 
of motion are described by Hamilton's principle: 

L{qit),q{t))dt^O. 

From the viewpoint of Dirac geometry, we start rather from the Hamilton-Pontryagin prin- 
ciple, which is given by the stationary condition of the action integral on the space of curves 
{q{t),v{t),p{t)), t e [ti,t2] in TQ®T*Q as 

L{q{t),v{t)) + {p{t)A{t) - v{t)) dt = 0, 



where q[t) = dq{t)/dt. The above variational principle can be restated by using the gener- 
alized energy El : TQ ® T*Q R as 

{p{t),m) - EL{q{t),v{t),p{t)) dt = 0. 

Keeping the endpoints q(to) and q{ti) of q{t) fixed, we can obtain the implicit Euler- 
Lagrange equations. 

The Lagrange-d'Alembert-Pontryagin Principle. Consider a mechanical system with 
kinematic constraints that are given by a constraint distribution Ag on Q. Then, the mo- 
tion of the mechanical system c : [^1,^2] — >• Q is said to be constrained if c{t) G Ag(c(t)) 
for all t, tl < t < t2. Assume that the distribution Aq is not involutive; that is, 
[X{q),Y{q)] ^ AQ(q) for any two vector fields X,Y onQ with values in Ag. Further, let L 
be a (possibly degenerate) Lagrangian on TQ and let F : TQ (B T*Q ^ T*Q he a.n exter- 
nal force field. The Lagrange-d'Alembert-Pontryagin principle for a curve {q{t),v{t),p{t)), 
tl < t < t2, in TQ <3) T*Q is represented by 

■ {p{t),q{t)) -EL{q{t),v(t),p{t))dt+ f ' {F{q{t),v{t),p{t)),6q{t)) dt 
tl Jti 

= 8 r L{qit),v{t)) + (p(t), q(t) - vit)) dt+ r {F{q{t), v{t),p[t)),5q{t)) dt 
Jti Jti 

= 

for a given variation Sq{t) E AQ{q{t)) and with the constraint v{t) e AQ{q{t)). Keeping the 
endpoints of q{t) fixed, it follows 

— ^p,Sqj + (^—-p,Svj + {Sp,q-v)dt + J F{q,v,p)Sqdt = (3.4) 

for a chosen variation dq{t) E AQ{q{t)), for all Sv{t) and Sp{t), and with v{t) E AQ{q{t)). 
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Proposition 3.1. The Lagrange-d'Alembert-Pontryagin principle gives the local expression 
of equations of motion for nonholonomic mechanical systems with external forces such that 

dL dL 
q = veAQ{q), p- — -F{q,v,p)eA°Q{q), (3.5) 

Proof. From equation (3.4), it reads 

(^ — -p + F{q,v,p),Sqj + -P^^vj + {Sp,q-v) dt = 0, 

which is satisfied for a given variation Sq{t) S AQ{q{t)), for all Sv{t) and Sp{t), and with 
the constraint v{t) g AQ{q{t)). Thus, we obtain equation (3.5). ■ 

Coordinate Representation. Let Q be a finite dimensional smooth manifold whose 
dimension is n. Choose local coordinates on Q so that Q is locally represented by an 
open set U C M". Then, TQ is locally given by {q\v') e C/ x M". Similarly T*Q may be 
locally denoted by {q^,Pi) G [/ x M". The constraint set Aq defines a subspace on each fiber 
of TQ, which we denote by AQ{q) C M" at each point q € U . If the dimension of AQ{q) 
is n — m, then we can choose a basis em+i{(i), 6^+2(5), • . • , en(<z) of AQ{q). Recall that the 
constraint sets can be also represented by the annihilator A°Q{q), which is spanned by m 
one-forms cj^, w^, . . . It follows that equation (3.5) can be represented, in coordinates, 
by employing the Lagrange multipliers /ia, a — \, ...,m, as follows: 

-w-FM.^^p)\^f 




where we employ the local expression w° ~ ujf dq^ 



4 Interconnection of Dirac Structures 

There exist various types of the power conserving interconnections between distinct dynam- 
ical systems such as two rigid bodies connected by ball-socket joints, mechanical systems 
connected by transmission gears, transducers between electric circuits and mechanical sys- 
tems such as D-C motors, fluid and rigid body interactions, etc. All the examples listed 
above may be represented by Dirac structures. In this section, we will show how distinct 
Dirac structures can be interconnected via an interconnection Dirac structure to yield an- 
other Dirac structure of a fully interconnected dynamical system. 

The Direct Sum of Dirac Structures. Before going into details on the case of n distinct 
Dirac structures, we will begin with the case of n = 2 to introduce the direct sum of Dirac 
structures on manifolds. Let ~I(M) be the set of (almost) Dirac structures over a manifold 
M.^ Let Ml and M2 be distinct smooth manifolds. To transition from systems on Mi 
and M2 to one on AI — Mi x M2, we shall define a direct sum associated to distinct Dirac 
structures. 

is called "Daleth", which is the fourth letter of the Semitic alphabet corresponding to "D". 
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Definition 4.1 (Direct Sum of Dirac Structures). Let Di e and D2 € ~I(M2) be 

given distinct Dirac structures on Mi and AI2 respectively. Define the Dirac sum of Di 
and D2 by 

Di ® D2{mi,m2) = {(K, t-a), T^, pr^, (ai) + T^, pr^^.^aa)) e T„M x r,;M| 

{vi,ai) e -Di(mi), {v2,a2) e L'2(™2)},Vm = (mi,TO2) e M 

We now show that ® is a mapping from ~I(Mi) x ~I(Af2) — > ~1(-A^)- 

Proposition 4.2. //Di e ~I(Afi),D2 e 1(^2), t/ien Di ® D2 e ~\{M). 

Proof. Let m — {mi,m2) & M and let (w,a) G Di ® D2{m). Notice that dim((Di © 
D2)im)) = dim{Di{mi)) + dim(L»2(™2)) = dim(Mi x M2) for m = (7711,7712) G M. By 
definition, t; = (wi,7'2) G TrnM and a = T,*^ pr^^j^ (ai) + prji/^ (02) for some ai G T,* ^Afi 
and a2 G T^^M2 such that (ui,q;i) G 151(7711) and (w2,Q:2) G D2{m2)- Then, one has, for 
each 771 G M, 

{{{v,a), {v,a))) = 2{a,v) 

= 2(T;„ prj,^^ (ai) + T,; pr^^^ (^2), v) 
= 2{ai,TmprM^{v)) + 2(a2, T™ pr^^^^ (w)) 
= 2(ai,?;i) + 2(02,-^2) 
= 0, 

since {vi,ai) G Z3i(777i) and (w2,a2) G -02(^2). Thus, noting that 771 is arbitrary, one can 
prove the theorem by Lemma 2.1. ■ 

The Dirac Sum of Induced Dirac Structures. Let Qi and Q2 be distinct configuration 
spaces. Let Ag^ C TQi and Aq^ C TQ2 be constraint distributions, and we can define the 
induced Dirac structures Di and D2, where we assume that Aq^ ^ Aq^ and Aqj^iHAq^ — 0. 

Let Q — Qi X Q2 he an extended configuration manifold and Aq = Ag^ x Ag^ C 
TQ — TQi X TQ2 be a distribution on Q. As before, define the induced distribution At-'q 
on T*Q by 

AT-g = (r7rg)-i(Ag)crr*g, 

where Tttq : TT*Q — > TQ is the tangent map of the cotangent bundle projection ttq : 
T*Q Q, while the annihilator of A^^g is defined by, for each {q,p) G T*Q, 

A^g(<Z,p) - {a G Tl^,pjT*Q \ {a,w) = for all w G AT.g(g,p)}. 

Proposition 4.3. Let be the canonical symplectic structures on T*Qi and il^ : TT*Qi — > 
T*T*Qi be the associated bundle maps. Let pr^-^g^ : T*Qi x T*Q2 T*Qt, i = 1, 2 be the 
natural projections'^ . Denote by Di the Dirac structures on T*Qi induced from Aq. C TQi. 
Then, the Dirac sum Di © D2 can be represented by the induced Dirac structure on T*Q 
from Aq = Agj x Ag^ C TQ, which is given by, for each {q,p) G T*Q, 

Di®D2{q,p) = {(u;,a) G T(,,p)r*g x Tl^^^-^T*Q \ w G AT.g(g,p) 

and a-{ni(B Vl2)\q,p) ■ w G A?,.g(g,p) }, (4.1) 

where {ili(Bi^2)^ '■ TT*Q — > T*T*Q is the associated bundle map with Qi(S^2 — pr^'gi ^1 + 
pr^.g^ ^2- 



*Note that we may identify the cotangent bundle T*Qi X T*Q2 with T*(Qi X Q2). 
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Proof. Let {w,a) G -Di © -D2(g,p) for each {q,p) e T*Q. Let prTT'Q* : TT*Q, ® 
T*T*Qi — > TT*Qi, i = 1,2. Now, we may decompose (w,a) as «; = (■U'i,W2) and 
a = ^(g,p)PrT*Qi(ai) + ^(g,p)PrT*Q2('^2) such that (wi.ai) € Di(gri,pi) and (■u;2,a2) S 
D2{qi-,Pi)- Define the distributions At*Qi = prry,Q^(Z)i) and At'Qs = 'i)Vj,j,,Q^{D'2). 
Then, It follows that ai — U\{qi,pi) ■ wi € A^,Q^{qi,pi) and a2 — Q2i'l2,P2) ■ wi € 
^T*Q2(927P2), where wi e AT-Qi(gi,Pi) and W2 € At*Q2(92,P2)- Noting that (Qi © 
f^2)^g,p) = © n\{q2,p2) and w = (^1,^2) € AT*Qi(gi,Pi) x AT*Q2(gi,P2), one 

has a - (r^i © 172)^(g,p) • w e A^^.g^ (91,^1) x A^^.g^ (g2,P2)- Additionally, Titq{1^t'q, x 
Ar-Qa) = Aqj x Aq2 = Aq C TQ. It follows on each fiber Tf^q^p)T*Q that At-Qi x At-Q2 = 
Ay.Q and it is obvious that A^.g = A^.g^ © A^.g^. Putting these all together gives us 
w e Ax'Qiqjp) and a - (fJi © rJ2)^(g,p) • w e A^.g(g,]3), for each {q,p) G T*Q. Since 
ill © f22 is the canonical symplectic form of T*Q we find Di © D2 is the induced Dirac 
structure from Ag. ■ 

An Interconnection Dirac Structure. In order to intertwine distinct Dirac structures, 
we will introduce a special Dirac structure called interconnecting Dirac structure. A simple 
case may be given by a distribution and its annihilator. 

Suppose there exists a constraint distribution Ac C TQ due to the interconnection, such 
that 

Ag / Ac and Aq n Ac ^ 0, 

where coq is a one-form on Q associated with the interconnection. The annihilator Ag C 
T*Q is defined by 

A°(<7) = {/ € T*Q \{f,v)=0\v€ Ac(g)}. 
We can define the lifted distribution Ac on T*Q by 

Aint = (r7rQ)-i(Ac)cTT*Q, 

where Tttq : TT*Q TQ is the tangent map of the cotangent bundle projection ttq : 
T*Q — > Q, while the annihilator of Ajnt is given by, for each {q,p) G T*Q, 

Al,{q,p) = {ae Tl^,p)T*Q I (a,w> = for all id e A,^,{q,p)). 

Then, the interconnection Dirac structure associated to the distribution Aint on T*Q 
can be defined by, for each {q,p) e T*Q, 

Dintiq,p) = { iw,a) e T^g,p)T*Q x %^p)T*Q \ w G Ai„t(9,p) 

and aeA°„t(q,p)}. (4.2) 

Example: Two Particles Moving with Contact. Consider two masses moving with 
contact in a vector space V = M^, whose velocities are given by (wi,W2) G V- Since the two 
particles are in contact and their velocities are common, it follows 

{vi,v2) e Ac c V, 

where Ac = {{v\,V2) \ vi = V2} is a constraint subspace of V. This constraint is enforced 

through the associated constraint forces F2) € V* at the contact point, where Fi is the 
force associated with the i-th particle. This implies Newton's third law and it follows 

(Fi,F2)e A°cy*, 

where A° = {(-Fi, -F2) I = — -^2} is the annihilator of Ac. Note that we can define a Dirac 
structure on V by Dy = Ac x A°. Then, it follows that (-yi, Fi), {v2, F2) G Dy- 
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Remark. We can consider a more general interconnection Dirac structure than stated in 
the above, such as those given by gyrators, transformers and so on, which are expressed by 
two-forms or Poisson structures. In this paper, we wih mainly focus on the simple case of 
the interconnection Dirac structure Di^t induced from a given distribution Ac 

Strategy for the Interconnection of Distinct Dirac Structures. Before going into 
details, let us outline our strategy for interconnecting Dirac structures. Let Di and D2 be 
two distinct Dirac structures on distinct manifolds Mi and M2. We will interconnect Di 
and Z?2 by introducing an interconnection Dirac structure Di^t on M = Mi x M2 , by which 
constraints due to the interconnection are imposed. To do this, we first make the direct 
sum of Dirac structures, namely Di D2, which is a single Dirac structure on M. Then, 
by introducing the bowtie product cxi of Dirac structures on M, we define an interconnected 
Dirac structure as {Di (BD2) i>i -Dint- The bowtie product will be shown to make (~I(M), cxi) 
a commutative category with identity element TM® {0}. Under the condition that the fiber 
dimension of pTrpj^i{Di) n ]irj^]^,j{D2) is constant, it will be shown that (Di D2) ix -Dint is 
indeed a Dirac structure on M . 

The Bowtie Product. Let P — M x M and let t : M ^ P be a submanifold given 
as M := {(m, n) E P \ m — n}. Given Dirac structures Da, Di, on M, define the direct 
product of Dirac structures Da and Di, by Da ® Dt on P. Recall from Courant [1990] that 
one may pull back the Dirac structure Da® Di, C TP ®T*P on P to M as 



which is a Lagrangian submanifold of TM ® T*M. Assume that i*{Da ® Db) is a smooth 
sub-bundle by either of the following conditions: 

• {Da X Db)n {TM ® T*P) has constant dimension, 

• {Da X Db) n ({0} ® TM°) has constant dimension. 

Corollary 4.4. The subhundle L*{Da ® Db) is a Dirac structure on M. 

Let <: : 7\/ — > M be the map m 1—^ (m, m) (this is a diffcomorphism between M and M) . 
Then d := l o <; : M ^ P given by m — ^ {m,m) is called the diagonal embedding. Given 
Dirac structures Da, Db € ~\{AI), we may define the Dirac structure Da ® -Df, € ~\{M x M). 
Noting i*{Da ® Db) G ~I(M), it follows that 



Definition 4.5. Let Da and Db be Dirac structures over M . Then, we define a tensor 
product called the bowtie product of Da and Db by 



t*{Da®Db) : 



{Da®Db)n{TM®T*P\M) 
{Da®Db)n{{0}(BTM°) 



d* {Da ® Db) = ^* [i*{Da ® Db)] e 1{M). 



Da^xiDb^ {{v, a) £ TM ® T*M \ 3/3 e T*M 

such that {v, a + (3) £ Da, {v, -(3) £ Db}, 



(4.3) 



or alternatively given by: 



Dat^Db=d*{Da®Db). 



(4.4) 



Theorem 4.6. Under the assumption that {Da(BDb)r\{{0}(BTM°) has constant dimension, 
Da [X Db is a Dirac structure on M . 
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Remark. In this paper, we focus on how to use the bowtie product ix to make a Dirac 
structure for an interconnected Lagrange-Dirac dynamical system. While the bowtie prod- 
uct was initially proposed in Yoshimura, Jacobs and Marsden [2010] for this purpose using 
formula (4.3), it was revealed that an equivalent tensor product (denoted by M) was also 
developed in the context of generalized complex geometry (see Gualticri [2007] and Alek- 
scev, Bursztyn and Mcinrcnken [2009]) using formula (4.4).^ We will make explicit the 
equivalence of cxi and Kl (i.e. (4.3) and (4.4)) in the Appendix. 

Remark. For standard Lagrange-Dirac mechanics, we usually employ M — T*Q, Da is 
the canonical Dirac structure on M, and Df, = A © A° for some distribution A C TM 
is an interconnection Dirac structure. Though we will not explore this in detail, one of 
the particularly interesting interconnection Dirac structures is the case in which Dt can be 
developed by Lagrangians that are linear in velocities of the form L = (a, u), where u € TQ 
and a is a one-form on Q causing magnetic terms to appear. 

In order to show interconnection of Dirac structures is associative, we will use a special 
restricted two-form J7a induced by a Dirac structure D on M with a distribution A = 
pryj,j(£>), where pttm ■ TM © T*M TM; {v,a) i-> v. This two-form wih be employed 
in the proofs for theorems 4.9 and 4.11. 

Lemma 4.7. Let D e 1{M) and set A = pry^f(D). On each fiber of T^M x T*M at 

X € M there exists a bilinear anti- symmetric map Q/\{x) : A(a;) x A{x) — > M defined as 



This two-form was initially introduced by Courant [1990] for linear Dirac structures (see 
also Dufour and Wade [2004]). We can easily generalize it to the case of general manifolds 
since is defined fiber-wise, and hence we have that D{x) is a linear Dirac structures on 
each fiber T^M x T*M for all x € M. Here, we present a proof for completeness. 

Proof. We desire to show that equation (4.5) defines a unique bi-linear anti-symmetric map 

First we prove linearity. Let {vi,ai), (^2,0:2) S D{x) and let Ci,C2 £ M. Since D{x) is a 
linear space at each x S M we have {ciVi + C2V2, c\a\ -f 020:2) S D(x). By definition of fl^, 
it follows, for all w E A(x), 



Thus, ilA is linear in the first argument. Second, in order to prove that fl^ is anti-symmetric, 
take any {vi,ai), {v2,a2) G D{x). Since D is isotropic, we find 



^a{x){vi,V2) — {oLi,V2) such that (wi,ai) G D{x). 



(4.5) 



Q.a{x){ciVi + C2V2, w) = (citti -I- 0202, w) 

= ci{ai,w) + C2{a2,w) 

= Ci^a{x){vi,w) + C2flA{x){v2,w). 



(((ui,ai), (W2,a2))) 



{ai,V2) + {a2,vi) = 



and we compute 



^a{x){vi,V2) = {ai,V2) 



{0:2, Vi) 



^^a{x){v2,vi). 



Thus, we have shown anti-symmetry from which bi-linearity follows. 



^We appreciate Henrique Bursztyn for pointing out this fact in private communication. 
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Definition 4.8. Given a Dirac structure D G ~I(M), we call the two-form JIa of Lemma 
4-7 the Dirac two-form. 

Recall from equation (2.1) that, for each x £ M, D{x) may be given by 

D{x) = {{v,a) £ T^M X T*M \ v £ A{x), and a{w) = nA{x){v,w) for all w £ A{x)}, 

which may be also stated by 

D{x) = {iv,a) £ T^M x T*M \ v £ A{x), and a~Vt{x) ■v£ A°{x)}, 

where fl^lx) : T^M — >■ T*M is the skew-symmetric bundle map that is a natural extension 
of the skew-symmetric map n^^ix) : A{x) C T^M A*{x) = T*M/A°{x) C T*M, which 
is defined by {i}^^{x){vx),Wx) — i^A{vx,Wx) on A{x). 

Theorem 4.9. Let Da and Di, £ ~\{M). Let Aq — pirpj^^j[Da), Af, — pirpj^j^Di,) and let ila 
and rib be the Dirac two-forms for Da and D^ respectively. If Aa n A^ has constant rank, 
then Da cxi Di, is a Dirac structure with distribution Wrmi^a ^ = Aa Ci Af, and with 
the Dirac two-form {fla + f^b)|A„nAb- 

Proof. Let {v,a) £ Da t>< Db{x) for each x £ M. Then, there exists (i £ T*M such that 
{v, a -\- /3) £ Da{x), {v, —/3) £ Db{x) for each x £ M . This is equivalent to 

n\{x) - v-a- P£ A°(x) and nl{x) ■ v + f5 £ Al{x), for each x £ M, 

where v £ Aa{x) and v £ Ab{x). This means (Jl^ + ill){x) ■ v - a £ A°(a:) + Ag(x) 
and i; e Aa n Ah (a;). But A°(a:) + A^(x) = (A^ n Ab)°{x). Therefore, upon setting 
= + and Ac = Aa n Ab, we can write ^^^{x) ■ v ~ a £ A°(a;) and v £ Ac{x); 
namely, {v,a) £ Dc{x), where Dc is a Dirac structure with Ac and ilc- Then, it follows 
that Da fxi Db C Dc- Equality follows from noting that both Da cxi Db{x) and Dc{x) are 
subspaces of TxM x T*M with the same dimension. ■ 

Corollary 4.10. If fib = 0, then it follows that Db — A;,® A^ and also that Dc = Da ix: Db 
is induced by the constraint Aa A;, and the two-form Via . 

Theorem 4.11. Let Da, Db, Dc £ ~\{M) with distributions A a = prj^j^j{Da), Ab = prrpj^j[Db) , 
and Ac — prrpj^j(Dc) . Assume that Aa H A^, Af, n Ac and Ac H A^ have constant ranks. 
Then the bowtie product Cxi is associative and commutative; namely we have 

{Da ix Db) cxsDc^ Dat^ (Db CXI Dc) 

and 

Dai=<l Db ^ DbCXI Da. 

Proof. First we prove commutativity. Recall that any Dirac structure may be constructed 
by its associated constraint distribution A = prp^(£') and the Dirac two-form IIa- Let 
fla, i^b,and Oc be the Dirac two- forms corresponding to Da, Db, and Dc respectively. Then 
we find by Theorem 4.9 that Da ex Db is defined by the distribution Aab — Aa H Ah and 
the Dirac two-form il^A^b — (^A„ + ^^Ab)|A„i,- By commutativity of -I- and H, we find the 
same distribution and the two-form for Db ex Da, we have Da cx Db = Db cx Da- 

Now we prove associativity. Let A(a&)& = WTMii^a tx Db) cx Dc) and Aa(bc) — 
prrpj^j[Da cx {Db cx Dc)) and it follows 



A(ah)c = (Aa n Ah) n Ac = Aa n (Ah n Ac) = Aa(hc)- 
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If ^A^ab)c ''^^'^ ^A„(i,c) are respectively the Dirac two-forms for {Da ixi Db) cx] Dc and Da t<\ 
{Df, ixi Dc), we find 

= (f^A„+f^A. +r!Aj|A„(,., 
= ^A„(,,). 

Tfius, we obtain 

{Da CXI -Db) txDc = Datx {Db >i Dc). 



Remark. We have shown that txi acts on pairs of Dirac structures with clean intersections 
to give a new Dirac structure and also that that it is an associative and commutative product. 
It is easy to verify from equation (4.3) that the Dirac structure I?e = TM © {0} satisfies 
the property of the identity element as D^t^ D — D — D \xi D^ ioT every D e ~\{M). Thus 
we find that (~I(M),ixi) is a commutative category (it lacks totality and an inverse element, 
so it is neither a monoid nor a groupoid). 

Interconnection of Two Distinct Induced Dirac Structures. Let Di and D2 be 

Dirac structures on distinct manifolds Qi and Q2 induced from constraint distributions 
Aqj C TQi and Aq^ C TQ2 as before. Let Ac be a given distribution on Q = QiX Q2 due 
to the interconnection of Di and D2. We have mathematical ingredients for interconnecting 
the induced Dirac structures Di and D2 through the interconnection Dirac structure Z3i„t 
on Q. Recall that an interconnection Dirac structure is given by Z?i„t = Ajnt © ^i'nf where 
Aint = {Tt^q)~'^{^c) C TT*Q is the constraint distribution associated to the interconnec- 
tion. RecaU that At-q^ = {Ttiq,)'^{1^q,) C TT*Qi and A2 = (TttqJ-HAqJ C TT*Q2. 
It follows from Theorem 4.9 that D — {Di D2) ixi Mnt may be constructed from the 
constraint A — (Ai x A2) n Ai^t and = i^i il2, where fii and rt2 are respectively the 
canonical two-forms on T*Qi and T*Q2- 

Proposition 4.12. Assuming A has constant rank, the interconnection of two distinct 
induced Dirac structures Di and D2 through Di^t = Ajnt ® A?^^, namely, 

D^{Di® D2) M Ant 

is represented by, 

D{q,p) = { {w,a) e T^g.p)T*Q x Tl^^p)T*Q \ 

w e A{q,p) and a — 51^(g,p) • w £ A°{q,p) }. 

This may be viewed as a Corollary to Theorem 4.9 and Proposition 4.1. 

Interconnection of n Distinct Dirac Structures. Let us consider the interconnection 
of n distinct Dirac structures Di, D2, ■ ■ ■ , D^ on distinct manifolds Mi, Af2, . . . , M„, in which 
n > 2. Recall that the Dirac sum ® is associative as {Da ® Db) (B Dc = Da ® (Db Dc), so 
that we may define the iterated Direct sum by 

N 

0A = A ©A© •••©£>„. 

i=l 
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By choosing the appropriate interconnection Dirac structure 

Ant e X • • • X M„) 

and assuming that rank(pry^,^(©Z?i) H pryj(j(Z?int)) is constant, we can simply develop the 
interconnected Dirac structure for distinct Dirac structures as 

^=|^0A^MA„t. (4.6) 

The condition that rank (prjn^j(0Z3i) H pr^^,^ (Ant)) is constant is equivalent to the as- 
sumption that rank(A) is constant, which can be shown by resorting to local coordinate 
expressions. 



A Link Between Composition and Interconnection of Dirac Structures. The 

notion of composition of Dirac structures was introduced in Cervcra, van der Schaft, and 
Banos [2007] in the context of port-Hamiltonian systems. Let Vi, V2 and Vg be vector spaces. 
Let Di be a linear Dirac structure on Vi ® Vs and D2 be a linear Dirac stucture on ® V2. 
The composition of Di and D2 is given by 

Di\\D2 = {{vi,V2,ai,a2)e{VixV2)(S{V{ XV2*) \ 

3{vs,as) eVs® V*, such that (wi, Ws, ai, a^) e Di, (-w^, W2, as, a^) G A}, 

where Vj* , V2 and V* denote the dual space of Vi , V2 and Vs ■ It was also shown that the set 
1)11102 is itself a Dirac structure on Vi x ¥2- In this section, we shall see the link between 
the composition and the interconnection of Dirac structures. 

The composition developed in Ccrvera, van der Schaft, and Banos [2007] is constructed 
on vector-spaces. So, we also focus on the vector space case in this paper. Let V = 
VixVsXVgX V2 and V — Vi x V2- Before going into details, we establish an important fact 
on the projection from V to V. Set : ^ — > to be the projection {vi,Vs,v's,V2) (wi, ^2). 
The mapping : V* V* dual to * is given by **(ai,a2) = ("1,0,0,02) G V* . 

Next, we recall the push-forward map associated to a Dirac structure. Let U and W be 
vector spaces. Let cp : U W he a linear map and let Djj be a linear Dirac structure on 
U. The push-forward of D{j to W hy (p is the set 

Dw = {Mv),a) eW®W* \veU, aeW\ {v,ip*{a)) G Du}- 

It is worth noting that the push- forward of a Dirac structure ip*Dij is itself a Dirac structure 
D\y; namely, we write it as 

ip*Du — Dw 

See, for instance, Bursztyn and Radko [2003] and Yoshimura and Marsdcn [2007b]. 

We are now ready to provide the link between the bowtie product and the composition 
of Dirac structures. 

Proposition 4.13. Lef Ajnt = {(wi, v^, -w^, ^2) e V'}- ^et Ant = Aint(SA°^^. * : V ^ 
V to be the projection (vi, Ws, Wg, ^2) '-^ {vi,V2). For linear Dirac structures Di on V\ x Vs 
and D2 on Vs XV2, set 

^{Di®D2)^ Ant 

and set also 

= AHA. 

Then, one has 

D\\ = ^-.I?^. 
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Proof. First, we observe that A?jj = {{0,as,as,0) G V*}. We also observe 

by definition of the push-forward map. Using the facts that ^{vi,Vs,v'g,V2) = {vi,V2) and 
**(ai,a2) = (ai,0,0,a2) e V*, 

={{vi,V2,ai,a2) \ ^Vs,v'^ e Vs such that {vi,Vs,v'^,V2,ai,0,0,a2) G D^}. 

Further, it follows from the definition of the bowtie product that 

^E'*£>lx ={{vi,V2,ai,a2) \ e Vs and 3(3 € V* such that 

(i'i,'(;s,i;^,i;2,ai + /3i , , /3g , q;2 +^2) e Z^i 0-02, 

{Vi,Vs,v's,V2,-Pl,-Ps,-l3s,-l32) G Ant}- 

Utilizing the fact that {vi,Vs, v'^, V2, — /3i, —(is, —P's^ —(^2) G Ant if and only if Vs = —v'^ and 
/3s = I3'g,(3i = 0, /32 = 0, we may restate the above as 

^*Dt^ ={{vi,V2,ai,a2) \ e Vs,as e V* such that 

{vi,Vs,-Vs,V2,ai,as,as,a2) G Di ®D2}. 

Finally, we have 

**-Dtxi ={{vi,V2,ai,a2) \ e Vs,as e such that 

(t;i,t;s,ai,as) e A, (-Ws,W2,Q!s,a!2) e A}, 
but this is nothing but the composition of Dirac structures, namely, D\\. ■ 

5 Interconnected Lagrange-Dirac Dynamical Systems 

Given distinct mechanical systems, one can develop dynamics of each mechanical system 
in the context of both the Lagrange-d'Alembert-Pontryagin variational structures and the 
induced Dirac structures over separate configuration manifolds. In this section, we will show 

how dynamics of an interconnected systems can be formulated in the context of the induced 
Dirac structures as well as in the context of the variational structures. 

Variational Structures for Interconnected Mechanical Systems. Let us introduce 
an interconnection of the Lagrange-d'Alembert-Pontryagin principles for distinct mechanical 
systems. 

Definition 5.1. Let Qi be distinct configuration manifolds corresponding to distinct me- 
chanical systems i = 1, ...,n. Let E^. : TQi ® T*(5i — >■ R be the generalized energies associ- 
ated to given Lagrangians (possibly degenerate) Lj on TQi and Fi : TQi © T*Qi T*Qi be 
external force fields. Let Aq. be a distribution on Qi for i = 1, n. 

The Lagrange-d'Alembert-Pontryagin principles for curves {qi{t),Vi{t),pi{t)), ti < 

t <t2 in TQi ® T*Qi is denoted by 



S {pr{t),qi{t)) - EL,{q^{t),Vi{t),p^{t))dt+ / {Fi{qi{t),Vi{t),p^{t)),^qi{t))dt 



= 5 Li{qi{t),Vi{t)) + {pi{t),qi{t)-Vi{t)) dt+ / {Fi{qi{t),Vi{t),Pi{t)),5qi{t)) dt 






= 



(5.1) 
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for Sqi{t) e AQ^{qi{t)) and with qi[t) e I^Q^[qi{t)). 

The interconnection of distinct Lagrange- d'Alembert-Pontryagin principles 

may he given by imposing some distribution Ac on the extended configuration manifold 
Q = Qi X ■ ■ ■ X Qn such that 

(gi,---,9n) e Ac(gi,...,g„) and © • • • © i^„(g„, i;„,p„) e A°((7i, g„). (5.2) 

Proposition 5.2. It follows from equations (5.1) and (5.2) that, by using variations Sqi{t) 
vanishing at the endpoints, one can obtain 

dL dL 
q^ = € Ag^iqi), p, - ~ Fi e A°Q.{q,), p^^^, i = l,...,n. (5.3) 

Proof. By direct computations and integration by parts, it follows that 



Li{qi,Vi) + {piAi - Vi) dt+ {F^, Sqi) dt 

+ Fi,5qi) + {5pi,qt - Vi) + {r>., - "1^-% ) dt 

= 

for 5qi € Aq. ((/i) and with qi G lS.Q^{qi). Imposing the interconnection condition (5.2), we 
obtain implicit Lagrange-d' Alembert equations for the interconnected system. ■ 

An Alternative Lagrange-d'Alembert-Pontryagin Principle. Here, we show that 
there exists an alternative Lagrange-d'Alembert-Pontryagin principle. Let Q ^ QiX - ■ - x Qn 
and let q = (gi,...,g„) G Q, {q,v) = (gi, (j„, wi, ?;„) G TQ = TQi x ■■■ x TQ„ and 
{q,p) = (gi, ...,(7„,pi, ...,p„) e T*Q 9^ T*Q^ x • • • x T*Q„. 

Define the generalized energy for the interconnected system by El = -Eli + • • • + i?L„ : 
TQ®T*Q — M, where L = Li + - ■ ■ + Ln ■ TQ — )• M is the Lagrangian for the interconnected 
system. Set the constraint distribution A = (Ag^ x • • • x Ag^) n Ac C TQ. 

The Lagrange-d'Alembert-Pontryagin principle for a curve {q{t),v{t),p{t)), ti < t < t2 
in TQ®T*Q may be given by 



(5.4) 

L{q{t),v(t)) + {p{t),q{t) - v{t)) dt = 0, 



{p{t)A{t))~EL{q{t),v{t),p{t))dt 
ti 
tl 

for the chosen variation 5q{t) G A{q{t)) C Tq(t)Q: for all 5v{t) and 5p{t), and with q{t) E 
A{q{t)) C r,(t)Q. 

Keeping the endpoints of q{t) fixed, it follows that one can simply derive the Lagrange- 
d'Alembert-Pontryagin equations as follows: 

A / N • dL . „ , - dL 
q = vEA{q), p- — eA°{q), P ^ j- ■ 
oq ov 

Proposition 5.3. The Lagrange-d'Alembert-Pontryagin principles in (5.1) for curves {qi(t), 
''^i{t) , Piit)) , tl < t < t2 in TQi © T*Qi together with the interconnection condition in (5.2) 
is equivalent with the Lagrange-dAlembert-Pontryagin principle in (5.4). 
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Proof. It follows from the Lagrange-d'Alembert-Pontryagin principles (5.1) that one can 
obtain equation (5.3) as 

dL dL 

4; e AQ,(gj), p^- - F, e A°Q^{q,), = i = l,...,n. 

Imposing the interconnection condition (gi , . . . , (7„ ) G Ac (gi ,...,(?„ ) and Fi (gi , , pi ) ® • • • ® 
e A°(gi, ...,(j„), one has 

(gi,...,(j„) = {vi,...,Vn) e A{qi,...,qn), \pi - -^,...,_p„ - j e A (gi,...,9„) 

and 



{Pl, ■■■,P2) 



\ dvi ' ' dv: 



where we note that the final distribution is given by A = (Ag^ x • • • x Aq^J n Ac C TQ 
and hence that its annihilator is denoted by 

A°(qi, g„) = A^^(gi) x • • • x A^Jg„) + A°(gi, g„). 



Theorem 5.4. T/ie following statements are equivalent: 

(i) A curve {q{t),v{t),p{t)), ti <t <t2 in TQ (BT*Q satisfies the Lagrange-d'Alembert- 
Pontryagin principle: 

6 [ ' {p{t)A{t)) - EL{q{t),vit),p{t)) dt 



^8 r L{q{t),v{t)) + {p{t),q{t) - v{t)) dt = 
Jti 

with respect to variations Sq (z A and Sv, 5p arbitrary and with fixed endpoints of q{t). 

(ii) Curves {qi(t) , Vi{t) , pi(t)) , ti < t < t2', i = l,...,n satisfy the Lagrange-d'Alembert- 
Pontryagin principles 

S f {p^^t),q^{t)) - EL{q^{t),V^{t),p,{t))dt+ f {F,{q,{t) , v,{t) , p,{t)) , 5q,{t))dt 

= 6 f L,{q,{t),v,{t)) + {p,{t),q,{t)-v,{t))dt+ f {F,{q,{t),v,{t),p,{t)),5q,{t))dt 
Jti Jti 

= 0, 

under variations 5qi{t) G AQ.[qi(t)), arbitrary Svi{t) and Spi(t), together with qi{t) S 
AQ.[qi(t)) and fixed endpoints of qi{t), such that the force vector F{q{t),v{t),p{t)) = 
Fii(li{t)iVi{t),pi{t)) e ••• © Fn{qn{t),Vn{t),Pn{t)) G and the velocity vector 

= (9i(0. -Anit)) e Tq(t)Q satisfy 

m e Ac(g(t)) and F{q{t),v{t),p{t)) e A°(g(i)), 

where Fi : TQi ®T*Qi — >■ T*Qi denote force fields. 



5 Interconnected Lagrange-Dirac Dynamical Systems 20 



Interconnected Lagrange-Dirac Dynamical Systems. Next, consider the intercon- 
nection of distinct Lagrange-Dirac dynamical systems corresponding to the above variational 
structures. To do this, we first consider the interconnection of distinct Dirac structures Di on 
T*Qi, i = 1, n, each of which is induced from given distinct distributions Ag^ C TQi. Let 
{E^., Di,Xi,Fi) hen distinct Lagrange-Dirac systems, where E^. denote the generalized en- 
ergies on TQi®T*Qi associated to given Lagrangians Li on TQi, : TQi®T*Qi — )■ TT*Qi 
are the partial vector fields, and the external force fields Fi : TQi ®T*Qi — > T*Qi. 

Recall that Fi induce the horizontal one-forms Fi on TQi T*Qi such that {Fi, Wi) = 
{Fi,TpTQ^(wi)) for all Wi G T{TQ, ® T*Q,), where prg^ : TQ, © T*Q, Qi are natural 
projections. 

As before, for an interconnected system, let Q = Qi x • • • x Qn and Ac be a given 
distribution on Q associated to the interconnection. Let A — (Aq^ x • • • x AQ^)n A^, C TQ. 
Then, define an interconnection Dirac structure by I?int = Ajnt © A?^^. as in (4.2), where 

Aint = (T7rQ)-i(Ac). Set a generalized energy El = El, ^ h El,, : TQ © T*Q M 

associated to the Lagrangian L = Li + ■ ■ ■ + Ln : TQ M. Furthermore, define a partial 
vector field on TQ © T*Q by X = Xi © ••• © X„ : TQ © T*Q TT*Q. Define also 
an external force field by F = Ti © • • • © : TQ ®T*Q T*Q, which induces the 
horizontal one- form F = Fi ® ■ ■ ■ ® F^ ow TQ ® T* Q such that {F, w) = {F, Tprgiw)) for 
ah w G T{TQ ®T*Q), where prg : TQ ® T*Q ^ Q is the natural projection. Recall from 
(4.6) that the interconnection of distinct Dirac structures Di,...,Dn through Di^t is given 
by L» = (L»i © • • • © D„) cx] A^^. 

Definition 5.5. The interconnected Lagrange-Dirac dynamical system is a quadru- 
ple {El, D, X) that satisfies 

{X,dEL\Tp)eD, 

where P = FL(A). Using local coordinates q — (gi,...,(j„) Cz Q, v = (wi,...,u„) G TQ, and 
p = {pi, ...,Pn) G T*Q, one has 

A / N • dL ^ « o . N dL 
q^veA{q), p- — -FeA°{q), P^g^- 

Proposition 5.6. The following statements are equivalent: 
(i) A curve {q{t),v{t),p{t)) G TQ®T*Q, ti<t<t2 satisfies 

(((7(i),p(t)),d£;i(g(i),«(i),p(<))|T(,(,),^(,,)P) G D{q{t),p{t)), 
where {q{t),p{t)) = ¥L{q{t),v{t)) G P. 
(a) Curves {qi{t),Vi{t),pi{t)), ti < t < t2 satisfy 

{{m,P^{t))AdEL^-F,){qmMt),P^mT,,^^,,.,^,,y,P.) ^ D,{q^{t),pm, 

where {qi,pi) G Pi = ¥Li{AQ.) , such that 

q{t) = {qi{t),...,q„{t)) G Ac(g(<)) and F{q{t),v{t),p{t)) G A°Q{q{t)), 

where F{q{t) ,v{t),p{t)) ^ Fi {qi (t) ,vi{t),pi{t)) ® ■ ■ ■ ® F„ (g„ (t) , w„ (t) , p„ {t) ) . 
Proof. Recall that D — {Di © • • • © _D„) txi D-mt, the condition 

(((Z(t),p(i)),di;(<7(t),«(i),p(<))|T„(„,^(,„p) G D{q{t),p{t)), for all ti < t < t^, 
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is equivalent to the existence of some covector 



such that 



and 



dL 

q,P,~-g^ - F,v + w ] e £'i(gi,pi) © • • • ® Ai(gn,Pn), (5.5) 



{q,P,F,~w) e Dint{q,p)- (5.6) 
It follows from the condition (5.5) that 

qi,p,,-- Fi,Vi + w.i] eDi[q^,pi), i = l,...,n, 



dqi 

and also from condition (5.6) it follows that q G Ac{q),w — and F G ^°{q)- Furthermore, 
noting dL/dqi = dLi/dqi, it follows that the curves {qi{t),Vi{t),pi{t)), ti < t < t2 satisfy 
the conditions 

({Ut),Mt))A<iE,-F,){q,{t),v^{t),p,{t))\ ) eD,{q,{t),p,{t)), 1 = 1,..., n, 

together with {qi{t),pi{t)) = VLi{qi{t),Vi{t)) e Pi. One can easily check the converse. ■ 

Thus, we obtain the following theorem: 

Theorem 5.7. The following statements are equivalent: 

(i) A curve {q,v,p)(t), ti < t < t2 in TQ (S T*Q satisfies 

{{q,p){t),<iE{q,v,p){t)\Tp) e Diqit),pit)), 

where {q{t),p{t)) = {q{t),{dL/dv){t)). 

(a) Curves {qi,Vi,pi){t), ii < t < ^2 satisfy 

((g,(t),p,(t)),(di?, -F0(9.(i),«.(0,P.(i))|T„^,„„^(,„pJ e D.,iq,it),p,it)), 

where {qi,Pi = {qi, dLi/dvi) and Fi is the horizontal one-forms on TQi T*Qi of 
force fields Fi : TQi (BT*Qi — > T*Qi, such that the condition of interconnection holds 

F{q,v,p){t) G A°((z(t)) and q{t) G A,((7(<)). 
(Hi) A curve {q,v,p)(t), ti < t < t2 satisfies the Hamilton-Pontryagin principle: 

i-t-2 

L{q{t),v{t)) + {p{t), q{t) - v{t))dt = 

ti 

with respect to variations 6q(t) G A{q{t)) and Sv, Sp arbitrary, with q{t) G A{q{t)) and 
with fixed endpoints of qit). 

(iv) Curves (qi,Vi,pi)(t), ti < t < t2 satisfy the Lagrange-dAlembert-Pontryagin princi- 
ples: 

6 L,{q,{t),v,{t)) + {p,{t),q,{t)~v,{t))dt+ {F,{t),Sq)dt = 0, 



where 5qi{t) G /S.Q^{qi{ty), and Svi{t) , Spi{t) are arbitrary, qi(t) G AQ.{qi{t)) and with 
fixed endpoints of q{t), such that the interconnection condition for the force vector 
F{q{t),v{t),p{t)) = Fi{qi{t),vi{t),pi{t)) ® ■ ■ ■ ® F„{qnit),Vnit),Pnit)) G r*(j)Q and 
the velocity vector q{t) = {qi{t), ...,g„(i)) G Tq^^^Q holds as 

q{t) G Ae(g(0) and F{q{t),v{t),p{t)) G A°(g(0). 
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6 Examples 

In this section we provide several examples of how interconnection can be carried out with 
the bowtie product. We have chosen rather simple examples, however the advantage of 
using the bowtie product is based on the procedure of tearing and interconnecting systems, 
which enables us to treat more complicated systems in a hierarchical structure. In future, 
we will seek more complex uses of the bowtie product. 

(I) Example: Mass-Spring Mechanical Systems. 

Let us consider an illustrative example of a mass-spring system in the context of the in- 
terconnection of Dirac structures and associated Lagrange-Dirac systems. In this example, 
there exists no external force. Let and ki be the i-th mass and spring (i — 1, 2, 3). 



Xi 



X2 



X3 



Figure 6.1: A Mass-Spring System 



Tearing and Interconnection. Inspired by the concept of tearing and interconnection 
originally developed by Kron [1963], the mechanical system can be torn apart into two sub- 
systems as in Fig 6.2, each of which can be regarded as a modular unit of the whole system. 
Let us call each disconnected system a "primitive systerri^ by following Kron [1963]. Note 
that the procedure of tearing inevitably yields interactive boundaries^ associated to the two 
disconnected subsystems, where it induces the following condition: 

/2+/2=0, ±2^X2. (6.1) 

Without the above condition, there exists no energetic interaction between the disconnected 
systems. In other words, the original mechanical system can be reconstructed by intercon- 



Sub-system 1 
Xi 



Tearing 



X2 



"'1 AAA/^ 



7712' 



Sub-system 2 
X2 X3 



f2 



Figure 6.2: Torn-apart Systems 



necting the subsystems through the continuity conditions. 

The continuity conditions in (6.1) imply the continuity of power flow; namely, the power 
invariance holds as 

V2+V2^ 0, 

where V2 — {f2,V2) and V2 = {f2,V2)- Needless to say, equation (6.1) can be understood as 
the defining condition for an interconnection Dirac structure. 

®Such an interactive boundary is called a "port" in circuit theory (see, for instance, Chua, Dcsoer and 
Kuh [1987]). 
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Primitive Lagrangian Systems. Let us consider how dynamics of the disconnected 
subsystems can be formulated in the context of Lagrange-Dirac dynamical systems. 

The configuration space of the subsystem 1 may be given by Qi = R x M with local 
coordinates {xi,X2), while the configuration space of the subsystem 2 is = x R with 
local coordinates {x2,X3). We can naturally define the canonical Dirac structures Di £ 
~\{T*Qi) and D2 G ~\{T*Q2) in this example. For Primitive System 1, the Lagrangian 
Li : TQi ^ M is given by 

Li{xi,X2,Vi,V2) = ^mivf + \^m2v\ - ^kixj - ^^(2:2 - Xi)'^, 
while the Lagrangian L2 '■ TQ2 — > M for Primitive System 2 is given by 

L2{X2,X3,V2,V3) = ^TOs^'s ~ \h{x3 " ^2)^- 

Then, we can define the generalized energy Ei on TQi(BT*Qi by Ei{xi, X2, vi^ V2, pi, P2) = 
PiVi + P2V2 — Li{xi,X2,vi,V2) and also define the generalized energy E2 on TQ2 ® T*Q2 

by E2{X2,X3,V2,V3,P2,P3) = P2V2 + P3V3 " ^2 (^2 , , ^2 , W3) ■ 

Though the original system has no external force, each disconnected primitive system 
has an interconnection constraint force at the interactive boundary. When viewing each 
system separately, the constraint force acts as an external force. This is because tearing 
always yields constraint forces at the boundaries associated with the disconnected primitive 
systems, as shown in Fig. 6.2 

The constraint force fields at the boundaries Fd : TQi © T*Qi T*T*Qi and 
Fc2 : TQ2®T*Q2 ^ T*T*Q2 are given by F,i = (xi, xa, wi, i;2,,Pi,P2, 0, -/a, 0, 0) an 
Fc2 = (X2, 2:3, z)2,W3,P2,P3,-/2, 0,0,0). Further, let : TQi ® T*Qi ^ TT*Qi be the 
partial vector field, which is defined at each point (pi = miVi,p2 = ni2V2) G Pi = VL{TQi) 
as Xi(xi,X2,-yi,W2,Pi,P2) = {xi,X2,Pi,P2,ii,i2,Pi,P2)- Similarly, let X2 : TQ2®T*Q2 -J> 
TT*Q2 be the partial vector field, which is defined at each point {p2 = 0,P3 = 'm3V3) G 
P2 = ¥L{TQ2) as ^2(^2, 2:3, "2, ^3,^2,^3) = (^2, 2:3,^2, Ps, ^2, ±3,^2,^3), together with the 
consistency condition p2 — 0. 

Primitive System 1 : We can formulate dynamics of Primitive System 1 in the context of 
the Lagrange-Dirac dynamical system {Ei, Di, Xi, Fi) as 

(Xi,di;i|TPi-Fci)e-Di, 

which may be given in coordinates by 

±1=1)1, X2=V2, Pl = -kiXi - k2ixi - X2), P2 = k2{xi - X2) + f2, 

and with pi = mivi and p2 = m2V2- 

Primitive System 2: Similarly we can also formulate dynamics of Primitive System 2 in 
the context of the Lagrange-Dirac dynamical system {E2, D2, X2, F2) as 

{X2, dE2\TP2 ^ Fc2) e -D2, 
which may be given in coordinates by 

X2=V2, X3^V3, P2 = k3ix3 - X2) + f2, P3 = "^3 (3^3 " ^2 ) , 

together with P2 = and p3 ~ ni3V3 as well as p2 — 0. 

Recall that each primitive system is physically independent of the other, which implies 
that there exists no energetic interaction between them. Next, let us see how the discon- 
nected primitive systems can be interconnected through a Dirac structure. 
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Interconnection of Distinct Dirac Structures. Let Q = (5iX(52=IRxIRxIRx]Rbe 
an extended configuration space with local coordinates x = {xi,X2,X2,X3). Recall that 
the Dirac sum of the Dirac structures is given by Di © D2 on T*Q. Now, the constraint 
distribution due to the interconnection is given by 

A,ix) = {v e T,Q \ {u;q{x),v) =0}, 

where ujq — dx2 — dx2 is a one- form on Q. On the other hand, the annihilator A° <Z T*Q 
is defined by 

A°(x) = {/ = (A, /2, h, h) e T;Q I (/, v) - and e Ac(a;)}. 

It follows from this codistribution that /2 = — /2, /i = and /a = 0. Hence, we obtain the 
conditions for the interconnection as in equation (6.1); namely, /2 + /2 = and V2 = V2- 
Recall that the distribution Ai„t on T*Q is obtained by Ai„t = {Tttq)'^{/^c) C TT*Q, 
and hence the Dirac structure Dint can be defined as in (4.2). Further, recall that the 
interconnected Dirac structure D on T*Q is given by 

Z?=(i?l©i?2)^Ant. 

Interconnected Lagrange-Dirac Systems. Now, let us see how two primitive Lagrange- 
Dirac systems, namely, {El-^, Di, Xi, Fi) and {El^, D2, X2, F2) can be interconnected to be 
a Lagrange-Dirac dynamical system. 

Define the Lagrangian L : TQ — )■ M for the interconnected system by L = Li + L2, 
and hence the generalized energy is given by El = El^ + ■ TQ (B T*Q R. Let 
A = (TQi X TQ2) n Aint. Set a partial vector field by X ^ Xi® X2 : TQ (S T*Q ^ TT*Q, 
which is defined at points {xi,X2,X2,X3,Vi,V2,V2,V3,pi = miVi,p2 = "^2^'2,P2 = 0,P3 = 
'msv^) € P — Fi(A), where the consistency condition P2 = holds. 

Finally, the interconnected Lagrange-Dirac system is given by {El, D, X), which 
satisfies, for each {x,v,p), 

{X{x,v,p),dELix,v,p)\Tp) G D{x,p), 
and with {x,p) e P = ¥L{A). 

The Lagrange-d'Alembert-Pontryagin Principle. Alternatively, we can develop the 

interconnected Lagrange-Dirac system by the Lagrange-d'Alembert-Pontryagin principle: 

S / Li{xi,X2,Vi,V2)+Pl{xi-Vi)+p2{x2-V2) 
J a 

+ L2{X2,X3,V3) +P2{X2 - V2) + PsixS " "3) dt, 

for all 6x2 — 6x2, for all Sv and 6p, and with V2 — V2. 

Thus, we can obtain the dynamics of the interconnected Lagrange-Dirac system. 
(II) Example: Electric Circuits 

Consider the electric circuit depicted in Figure 6.3, where R denotes a resistor, L an inductor, 
and C a capacitor. 

As in Figure 6.4, we decompose the circuit into the disconnected two circuits, each of 
which we shall call "Primitive Circuit 1" and "Primitive Circuit 2", in which and ^2 
denote external ports. In order to reconstruct the original circuit in Figure 6.3, the external 
ports may be connected by equating currents across each. 
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Figure 6.3: R-L-C Circuit 
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Figure 6.4: Disconnected Circuits 



Primitive Circuit 1. The configuration manifold for circuit 1 is denoted by Qi = M.^ with 
local coordinates qi = (qj^, q^, qsi), where qn, qL and qsi are the charge fluxes through the 
resistor R, inductor L and port Si respectively. There exists the KCL constraint distribution 
C TQi in Circuit 1, which is given by for each qi = (qji, q^, (75J G Qi, 

AQi(qi) = {vi = {vr,vl,vsJ e Tq^Ql \vb.~vl + Vg^ = 0}, 

where vi — {vr,vl,vs-^) denotes the current vector at each qi, while the KVL constraint is 
given by its annihilator Aq^ . Then, we can naturally define the induced Dirac structure Di 
on T*Qi from Ag^ as before. 

The Lagrangian for Circuit 1, namely, Ci on TQi is given by 

^i{qi,vi) = ^-^i^'i' 

which is degenerate. Define the generalized energy on TQi®T*Qi by i?£i(gi, = 
{PiiVi) — Ci[qi,vi) on TQi ®T*Qi. Circuit 1 also has the external force field due to the 
resistor and the external port, Fi : TQi T*Qi — > M as 

Fi{qi,vi,pi) = {-RvB.)dqR + fsidqs,- 

So, we can set up the Lagrange-Dirac dynamical system {Eci, Di, Xi, Fi). 

Primitive Circuit 2. The configuration manifold for Circuit 2 is = with local 
coordinates (72 = (9S2i9C')j where qs2 is the charge flux through the port 52 and qc is the 
charge stored in the capacitor. The KCL constraint distribution is given by, for each (72, 

A2(92) = {V2 = {vs2,vc) G Tq2Q2 \ Vc - -^52 = 0} 

and so the KVL space is given by its annihilator A2(g2)- This gives us the Dirac structure 
D2 on T* Q2 ■ Set the Lagrangian £2 : TQ2 M for Circuit 2 as 

1 2 
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and set the generalized energy i?£2(g2, ^2,^2) = (P2, ^2) ->C2(g2, 'y2) onTQ2®T*Q2. Circuit 
2 has the external force field due to the resistor and the external port S2, namely, F2 : 
TQ2®T*Q2 as 

F2{q2,V2,P2) = -fs2dqs2- 
Then, we can formulate the Lagrange-Dirac dynamical system {Ec2, D2, X2, F2). 

The Interconnection Dirac Structure Set Q ~ Qi x Q2- The interconnection con- 
straint is given by 

Ac = {u = {vi,V2) e TQ I t>Si = VS2} 
By the tangent lift Ajnt = Tttq (A12), we can define the interconnection Dirac structure as 

Ant = Aint © A°„t, 

which is denoted, in local coordinates, by 

Ant = {('ji,g2,Pi,P2,ai,a2,u'i,W2) G TT*Q®T*T*Q \ 

QSi = qs2 = 0, wi = 0, W2 =0, ai + a2 G span(dgsi - rf^Sa)}- 

The Interconnected System. Now, we can define the interconnected Dirac structure 

by 

i? = pi©A)i>aAnt, 
which may be viewed as the induced Dirac structure from the constraint space 

A== (Ai X A2)nAc, 

which is given, in coordinates (gi, 92) = (<?fi, qL.qSt , qs2,qc), by 

A(gi,(?2) = {(wi, W2) = {vr,vl,vsi,vs2,vc) \ vr - vl + vc = 0, vs\ = vs2, vs^ = vc}- 

Set the Lagrangian for the interconnected system as £ = £1 + £2 and the external force 
field F = Fi ® F2 : TQ ® T*Q T*Q. Set also Ec = Ec, + Ec2- Let X ^ Xi ® X2 : 
TQ®T*Q — )■ TT*Q be a partial vector field. The interconnected Lagrange-Dirac dynamical 
system is given by the quadruple {Ec, D, X, F), which satisfies 

{X, {dEc ~ F)\tp) e D, 

where P = FL(A). 

(Ill) Example: A Rolling Ball on Rotating Tables 

Consider the mechanical system depicted in Fig. 6.5, where there are rotating tables via a 
gear and a ball is rolling on one of the tables without slipping. We assume that there is no 
loss of energy and non-slipping in the gear and also that the external torque is constant. Let 
li be moments of inertia for the table i = 1,2. Let us decompose the system into distinct 
three systems; (1) a rotating (small) table, (2) a rotating (large) table, and (3) a rolling ball 
on a table. 

Primitive System 1. For Primitive System 1 the configuration manifold is the circle, 
Qi = S^. The Lagrangian is the rotational energy of the system. 

ii(si,si) = -^sj, 

and set the generalized energy El-^ on TQi®T*Qi. We employ the canonical Dirac structure 
on T*Qi as 

Di = {(si,_Psi,asi, w^J I si = Ws^,ps^ + as, = 0} 
There is a constant torque acting on Primitive System 1, given by a map Fi : TQi®T*Qi — >■ 
T*Qi defined as i^(si, WsuPsi) = tcLsi, where t is constant in M. 
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Figure 6.5: A Rolling Ball on Rotating Tables without Slipping 

Primitive System 2. For Primitive System 2, the configuration manifold is also the 
circle, Q2 — and the Lagrangian is again the rotational energy 

L2{S2,S2) = 

and set the generalized energy El^ on TQ2 ® T*Q2- Again, we have the canonical Dirac 
structure 

D2 = {{s2,Ps2,(^S2,Ws.,) I S2 = Ws^,Ps^ +as^, = 0}. 

Primitive System 3. Primitive System 3 is a rolling sphere of uniform density and 
radius l.The sphere moves in space by changing its position and orientation relative to a 
reference configuration. The configuration manifold is given by the special Euclidean group 
Q3 = SE(3), which we parametrize as {R,u) where R G S0(3),m S M^. Following Marsdcn 
and Ratiu [1999], let /? be the set of points of the sphere in the reference configuration. For 
configuration {R, u) G Q^, a. point x G P is transformed into by the action {R, u) ■ x = 
(R ■ x) + u. The Lagrangian is given by the kinetic energy as 



L^(R, u, i?, u) — 




where ||i?a; + M|p = x'^ R^ Rx + 2x^ W u + v? . We use body coordinates such that the center 
of the sphere in the reference configuration is at the origin so that xdx = 0. Substituting 
these relations, the above Lagrangian is to be 

J (^x^ifkx + ii?^ dx. 

Setting ma = pdx — ^np and noting that XiXjdx = when i 7^ j, one finally obtains 

L3 = !^(tr(7?^7?)+^.2) 

and with the generalized energy El^ on TQ3 ®T*Qy,. 

Since the motion along the z-direction is constrained so that the ball does not leave the 
plane of table 2, we get the holonomic constraint 

A3 = {(i?,i*) |«3 = 0}. 

This yields the induced Dirac structure 

-D3 = {{5R, 5u, 6pR, Spu, a_R, a„, wr, Wu) e TT*Qs, ® T*T*Q3 \ 

Sus = 0,Su = Wu, SR = wr, SpR + aR = 0, Spu + Uu — Xdz for some A G M}. 
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Interconnection Constraints. Let Q = Qi x Q2 x Q^. In order to interconnection the 
three primitive systems, we need to impose the constraints due to the non slip conditions. 
First, consider the mapping — s~^Vs for Vg € T^S^. The map (f> sends the tangent 

fiber TgS^ to the Lie-algebra s — T^S^ « R, where e is the identity element on viewed 
as a group. From this point on, we do not explicitly invoke $ but use it to interpret TS^ 
as X M. The interconnection constraint between System 1 and System 2 is given by 

Ac,i = {(si, S2, i?, u)eTQ\si + S2 = 0} 

and with its annihilator 

^c,i — span(a;i) 

where Wi = dsi — ds2- This constraint ensures that the gears rotate (without slipping) at 
the same speed in opposite directions. 

Next, we consider the interconnection constraint between Systems 2 and 3. Note that 
the velocity of a point located at the bottom of the sphere is given by 

= RR^ 

Note also that a point rotating on the gear of system 2 with the axle taken to be the origin 
has velocity 

-S2 A fx 

S2 J'[y 
So the non slip condition between System 2 and 3 is given by 

Ac,2 = {(si, S2, R,u)eTQ\ i • {-RR^ • k + u) = -53 • wj, j • {-RR^ • k + w) = §2 • Ui}, 
where i,j,k are the basis on M.'^. Setting the interconnection constraint distribution 

A, = n Ae,2, 

let Aint ■— iTTrQ)^^{Ac) and with its annihilator A?^^. Then, the interconnection Dirac 
structure is given by 

Ant = Aint © A°nt. 

The Interconnected Lagrange-Dirac System. The Dirac structure for the intercon- 
nected system is given by 

D = iDi(BD2(S D3) M Ant 

and we use the Lagrangian L = Li + L2 + L3. The dynamics of the system are implied by 
the constraint: 

{X,dEL\Tp) e D, 

where X : TT*Q®T*T*Q TT*Q is a partial vector field, and El ■.TQ®T*Q ^Ris the 
generalized energy of L. This gives us the interconnected Lagrangian system {El,D,X). 
By Theorem 4.9, these dynamics use the constraint distribution 

A = (rQi®rg2® A3)nA„ 

and the canonical two- form on T*Q restricted to A as the Dirac two-form. Noting the 
annihilator is given by 

A° = A°-HA°i + A°2, 

we finally obtain the dynamics as 

„ dL 
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7 Conclusions and Future Works 

In this paper, we provided a notion of interconnection for Dirac structures and their as- 
sociated Lagrange-Dirac dynamical systems. The principal idea lies in the fact that the 
interconnection of distinct Dirac structures itself is a Dirac structure. We also clarified 
the interconnection of the associated Lagrange-Dirac dynamical systems with the varia- 
tional structure by the Lagrange-d'Alembert-Pontryagin principle. To do this, we developed 
the bowtie product, [xi, which serves us well in the case of interconnection by constraints. 
We then showed how the bowtie product can effectively carry the Lagrange-d'Alembert- 
Pontryagin principle to the interconnected system. Lastly, we demonstrated our theory 
with some illustrative examples of a mass-spring mechanical system, an interconnection in 
circuits, and a non-holonomic mechanical system. 

We raise the following topics as arena's for further investigation: 

• Applications to complicated systems such as guiding central motion problems, multi- 
body systems, fluid-rigid body interactions, passivity controlled interconnected sys- 
tems, etc. (see, for example, Littlcjohn [1983]; Feathcrstonc [1987]; Yoshimura [1995]; 
van der Schaft [1996] and Ortega, van der Scliaft, Maschke and Escobar [2002]). 

• Reduction and symmetry for interconnected Lagrange-Dirac systems (Yoshimura and 
Marsden [2007b], Yoshimura and Marsden [2009]). 

• The use of more general interconnection Dirac structures with nontrivial Dirac two- 
forms such as the one given with symplectic gyrators (in this paper, we restrict Dint 
to the form Ajnt © ^hit)- examples of such general interconnections in physical 
systems, see for instance Wyatt and Cliua [1977]; Yoshimura [1995] and Maschke, van 
der Schaft and Breedveld [1995]. 

• The integrability condition for the bowtie product of Dirac structures. As to the 
integrability condition for Dirac structures, see Dorfman [1993] and Dalsmo and van 
der Schaft [1998]. The link with symplectic categories (see Weinstcin [2009]). 

• Discrete version of interconnected Lagrange-Dirac mechanics; namely, discretizing the 
Hamilton-Pontryagin principle one arrives at a discrete mechanics version of Dirac 
structures (see Bou-Rabee and Marsden [2009] and Lcok and Ohsawa [2010]). 

A Appendix: The Bowtie Product tx] and the Dirac 
Tensor Product lEI 

We will show the equivalence between the bowtie product cxj and the Dirac tensor product 
M developed by Gualtieri [2007], where the definition of the tensor product Kl is given as 
follows. 

Definition A.l. Let Da,Di, E ~\{P)- Let d : P ^ P x P be the diagonal embedding in 
P X P. The tensor product of Da and Df, is defined as 

K 

where K ^ {(0, 0)}©{(-C, ] C e T* P} C T{P x P)®T*{P x P) and is the orthogonal 
complement in T{P x P) Q) T* {P x P). 
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Proposition A. 2. The bowtie product tx and the tensor product M of Da and are 
equivalent: 

where the notation = implies the fact that these structures are isomorphic. 
Proof. First note that 

K-^ = {{v,v)}®T*{P X P) 
so that K + K-^ = {{v,v)} ®T*{P X P). Then we see that 

{Da®Db)n{K+K^) = {{v,v,aA,aB) € T{P x P)®T*{P x P) \ {v,v,aA,aB) € Da®Db} 
and the equivalence class Da ^ Db is 

Da^Db^ {{v, V, a, 0) I 3/3 such that {v, v,a + P, eDa® Db}, 
where we chose the element {v, v, a, 0) as the representative for the equivalence class 

[{v, V, a, a')] = {{v,v,a + ^,a-OMeT*P} eDa^ Db- 
Splitting the expression Da Db into Da and Db. we obtain 

Da^Db = {{v, V, a, 0) e T{P x P)®T* {P x P) \ 3(3 such that {v, a+(3) e Da, {v, -(3) e Db} 
we see that the assignment ^ : {v, v, a, 0) i-)- {v, a) is a bijection and also that 

(f>{Da K Db) = {{v, a) e TP e T*P \ 3/3 such that {v, a + /3) e Da, {v, -/3) e Db} 
this is just the definition of Da xi Db- ■ 
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